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Abstract 

We establish a new model, which takes into account a dynamic (inertial) self-interaction 
of gravitating systems. The model is formulated by introduction of a new function depending 
on the square of the covariant derivative of the velocity four-vector of the system as a whole 
into the Lagrangian. This term is meant for description of both self-action of the system 
irregularly moving in the gravitational field, and back-reaction of the motion irregularities on 
the gravity field. We discuss one example of exact solution to the extended master equations in 
the framework of cosmological model of the FLRW type with vanishing cosmological constant. 
It is shown that accelerated expansion of the Universe can be driven by traditional matter 
with positive pressure (e.g., dust, ultrarelativistic fluid) due to the back-reaction of the gravity 
field induced by irregular motion of the system as a whole; this back-reaction is shown to be 
characterized by the negative effective pressure. 
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1 Introduction 



The observational fact that the Universe expands with acceleration is till now a puzzle for theo- 
retical cosmology. One can formulate three main ideas, on which the explanations of this phenomenon 
could be based. The first idea is that there exists some exotic substratum, Dark Energy, possessing a 
negative effective pressure (see, e.g., [HO [3] for review and references). The second one is focused on 
modifications of the so-called geometric sector of relativistic theory of gravity: one can mention, for 
instance, f(R) theory, the Gauss-Bonnet model, etc. (see, e.g., [U El [6] for review and references). 
The Lagrangians of such theories contain the invariants constructed with nonlinear combinations 
of the Ricci scalar and the Riemann and Ricci tensors. The third idea introduces interactions of 
new types between gravity on the one hand, and fields and matter on the other hand (see, e.g., 
[H 13 13 [10]). The Non- minimal Field Theory is the most elaborated trend in this direction, and the 
corresponding models can be described by introduction of cross-invariants into the Lagrangian, which 
contain all admissible convolutions of the Riemann, Ricci tensors and Ricci scalar in combinations 
with the field strength tensors (see, e.g., [8| 19} fTTJ) . 

The theory of dynamic self-interaction of gravitating systems, which we establish here, is in line 
with the third idea. The motivation of the dynamic extention of the gravitational theory has two 
aspects: mathematical and physical ones. Started from the mathematical point of view, one can see 
that pure geometrical objects, which we use in the Lagrangian of the gravity field (the Riemann, 
Ricci tensors and Ricci scalar), contain second order partial derivatives of the metric. Since the 
covariant derivative of the metric itself is considered to be equal to zero, there are no geometric 
invariants containing the first derivatives of the metric only. Moreover, the covariant derivative of 
the scalar field reduces to the partial derivative, the field strength tensors for the electromagnetic 
and gauge fields, being the skew-symmetric quantities, in fact do not contain the Christoffel symbols, 
the four- divergence of the vector potential is assumed to be vanishing due to the Lorentz gauge. 
Thus, the Lagrangians of the Einstein-Maxwell, Einstein- Yang-Mills-Higgs, etc. theories do not 
contain invariants, into which only the first derivatives of the metric enter. The situation changes 
essentially, when it is acceptable to introduce into the Lagrangian the invariants containing the 
covariant derivative of the velocity four-vector, attributed to the macroscopic motion of the system 
as a whole. In this case one deals with the non-vanishing first derivatives of the metric, and such 
a gravitational theory can be regarded as a modification of the Einstein theory of a new type. The 
mathematical arguments supporting such a modification can be supplemented by the physical ones. 
Indeed, an accelerated point-like electrically charged particle produces electromagnetic radiation and 
is influenced by the back-reaction force, proportional to the second derivative of its velocity |12j . 
This is a clear example of dynamic self-interaction fulfilled by means of the electromagnetic field, 
the particle acceleration being the necessary condition of such self-interaction. It is well-known 
also, that irregularities of a medium motion lead to the dynamo-optical phenomena, which can be 
described by including covariant derivatives of the velocity four-vector into the Lagrangian of the 
moving electrodynamic systems in combination with the Maxwell tensor |12j . The back- reaction 
of irregularly moving electrodynamic system on the gravity field was discussed in [13J. A question 
arises: is it possible to describe a self-acceleration / deceleration of a material system as a result of a 
gravitational back-reaction on its irregular motion? We show below that it is indeed possible, when 
the Lagrangian of a matter contains some complementary invariant constructed as a square of the 
covariant derivative of the velocity four-vector of a system as a whole. As a cosmological consequence 
of this theory we obtain exact solution, which describes accelerated expansion of the Universe filled 
with traditional matter with positive pressure in the case, when the lambda-term is vanishing. 
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2 Dynamic extention of the gravitational theory 



Let us consider the variation of the action functional 

R 



S[gik] = J d 4 x^g 



2k 



+ L (m) + F(^ 2 



(1) 



with respect to metric gi k . The determinant of the metric, g = det(gi k ), the Ricci scalar, R, and 
the Lagrangian of a matter, Lr m ), are the standard elements of the action in the Einstein theory of 
gravity [TJ]. The new element in ([1]) is the term F(\I/ 2 ), which is presented by a generic function of 
the argument \I/ 2 



ik 



m ik = WiU k 



(2) 



introducing a covariant derivative VjC4 of the time-like four-vector field U k . We assume that U k 
is normalized by unity, i.e., gi k U l U k = 1, and is interpreted as a velocity four- vector of the matter. 
More precisely, we consider U k as a time-like eigen-vector of the stress-energy tensor of the matter 
T/ fc m) (the definition given by Landau and Lifshitz) 



rp(m) __ 

1 ik — 



2 5 V~9 ^(m) 



-g 



Sg 



ik 



(3) 



The eigen-value W relates to the scalar of the energy density of the matter, and the stress-energy 
tensor of a matter can be now decomposed as 



T, 



(m) 



ik 



WUiU k + v.. 



ik 



V, 



ik 



PA ik + n 



ik ■ 



(4) 



Here Vi k is an anisotropic pressure tensor containing a sum of the Pascal (isotropic) pressure, P, 
and a non-equilibrium pressure tensor, LT^. The tensor = gik~UiU k is a projector. Thus, the 
first (algebraic) relation in (j3J) allows us to use one scalar W, one four-vector U k and one tensor 
Vi k instead of ten components . The alternative approach based on the Eckart definition of the 
macroscopic velocity four- vector involves an additional particle number four- vector N k , a heat flux 
four-vector q k and the corresponding thermodynamic relations, thus sophisticating the model. 

The tensor ^>ik can also be decomposed into the sum of irreducible parts: acceleration four- vector 
DU l , shear tensor cx^, rotation tensor u^, and expansion scalar G 



^ik = ViU k 
These basic quantities are defined as 
DU k = U m V m Uk 



UiDU k + cr ife + u ik + -A ik Q . 



a lk = -A™A£ (V m U n + V n U m ) - -A ik Q , 



(5) 



u lk ee -A?A n k (V m U n - V n U m ) , e 



V m C/ r ' 



The quantity \l/ 2 , the argument of the function F(\l/ 2 ), has the form 



* 5 



DU k DU k + a ik a ik + uj lk uf k + ^0 2 . 



(6) 



(7) 
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The variation of the action functional ([I]) with respect to the metric yields 

Rik — -R 9ik = « (r^f 1 + T ik j ■ 

The tensor is defined by (JHJ) and (Tj0). The new term T ik , defined as 

=0 F(M/ 2 )] 



2 5 



ik 



Sg 



ik 



(8) 



(9) 



can be represented in the following form 



% k = g lk F{^ 2 ) + 2 



d 2 F 



2^2 



U(i^k)m — ^(ik)U m 



+ 2 



(^ M ^ m + ^[mk^r) + u {l v k) e + u {l R k)m u m — (d + e)# 



(10) 



As usual, the symbols \l/(j m ) and ^[i m ] denote symmetrization and skew-symmetrization, respectively. 
Since the velocity four-vector is normalized, i.e., gikU l U k = 1, the variation procedure with respect 
to metric is supplemented by the variation of the velocity: 



sir 



- A 5 g ™ (y v 5\ + u q 5i 



(see, e.g., [15] for details). The sum of two tensors Tj> k and % k should be divergence-free due to the 
Bianchi identities 



V fe 



T, 



(m) 



0. 



(12) 



ik T" 2 ik 

To make the model complete, one should formulate constitutive equations for P and Ui k . Here we 
restrict ourselves by the ansatz, that the Pascal pressure is given by a barotropic function P = P(W), 
and the non-equilibrium pressure is a function of the energy density, acceleration four-vector, shear 
tensor and expansion scalar, i.e., Ui k =Ui k (W, DU l , a mn , 9). Thus, we obtain a self-consistent model. 



3 Application to the space-time of the FLRW type 
3.1 Key equations 

Let us assume the metric to have the following form 

ds 2 = dt 2 - a 2 (t) (dx 2 + dy 2 + dz 2 ) , (13) 

and the velocity four- vector to have only one component U k = 5*. We deal with the so-called 
synchronous frame of references. The tensor ^i k for this case 

%k = ViU k = -T° lk = -g ik = (g ik - UiU k ) , (14) 

is symmetric. The dot denotes the ordinary derivative with respect to time. The decomposition (jSJ) 
has now only one non- vanishing irreducible component 

m ik = \A ik Q = -A lk = H{t)A ik , (15) 
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and the scalar \1/ 2 reduces to 3 if 2 . The stress-energy tensor (jTOl) converts into 



2 ^°o°fc 



SI) H 



11 

dH 2 



(16) 



where a new unknown function f(H) = F(3H 2 ) is introduced. Surprisingly, the four- divergence of 
the tensor ffTBT) vanishes for arbitrary f(H), i.e., V fc ^[if] = 0. Thus, we obtain from (j!2p that 



/A: 



also, and the energy and momentum of the matter conserve separately. Clearly, the 



velocity four-vector U l = S l Q is the time-like eigenvector of this tensor, related to the eigen-value 

df 



^H) = f(H)-H^. 



(17) 



Consequently, for the model under consideration the velocity four-vector U l happens to be an eigen- 
vector for the total stress-energy tensor T-^ + T ik} the corresponding eigen-value being the sum 
W + uj{H). Other three eigen- values coincide and are equal to 



7T 



(1) 



7T 



(2) 



7T (3 ) = n(H) 



V ' 3 dH 2 



One can interpret the quantities w(H) and 7r(H) as an effective (inertial) energy density and effective 
(inertial) pressure, respectively, which are produced by irregularities of the macroscopic motion of 
the system as a whole. The signs of these quantities can be positive or negative depending on the 
choice of the function f(H). Let us note that if f(H) is linear, say, f(H) = 3AH, then H-^L = f(H), 

= 0, and this function disappears from the key equation, since u(H) = n(H) = 0. Clearly, for 
such a case the corresponding terms in the action functional 



J d^xy/^AQ = J d 4 xAd k (y^gU 1 



(19) 



gives a complete divergence and disappears at the variation procedure. 

The stress-energy tensor of matter (J3J) inherits the symmetry of the FLRW space-time, thus, 
Vik = —V(t)Aik, where V{t) can generally include both Pascal (P) and non-equilibrium (IT) parts, 
V(t) = P(t) + n(£). One can reduce the equations for the gravity field to the following two equations 



3H< 



K 



(20) 



(21) 



Using a differentiation of the first Einstein equation ( 1201) and excluding H from ( f2Tj) . we can obtain 
the standard balance equation for the energy and momentum of the matter 



W + 3H(W + V) = 



:22) 



which, indeed, does not contain information about the function f(H). 

Let us consider an exactly solvable model in which the lambda-term is considered to be vanishing 
and the equation of state is presented by the standard barotropic formula 



V= ( T -1)W. 



(23) 
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thus assuming that the Pascal (equilibrium) pressure, P, is proportional to the energy density W, 
and the non-equilibrium pressure IT is absent. As usual, for such an equation of state one obtains 
from fl22j 



Pit) 



(7 



W(t) = W 



a(t) 
a(t Q ) 



-37 



In order to obtain H(t) we have to solve the key equation 

3H 2 + K[Hf(H)-f(H)] = KW 
The acceleration parameter —q(t) 



a(t ) 



1-37 



(24) 



(25) 



aH 2 



1 + 



H 

IP 



can be presented in terms of f(H) and its derivatives as 

-Q[f(H(t))} 



1-7 



{3H 2 + k 


Hf'(H)-f(H)]} 


IE 2 


1 + f f"{H) 





(26) 



(27) 



This formula is a direct consequence of ( 12~TT) . ( 1201) and (1231) . When /(.ff) = 0, the acceleration 
parameter — q = 1 — ^ is negative for the traditional matter with 1 < 7 < |, as it should be, 
nevertheless, the situation changes significantly, when f(H) 7^ 0. The form of the function f(H) is 
a subject of modeling. In the next papers we hope to consider this function in context of fitting of 
observational data, but here we restrict ourselves by the analysis of one model only, which allows us 
to find an example of exact solution and to study the acceleration parameter explicitly. 



3.2 Example of exact solution 

We consider the function F(\I/ 2 ) to be quadratic function of its argument, i.e., the function f(H) = 
F(3H 2 ) has the form 

Kf(H) = a + f3H 2 + aH 4 . (28) 

According to (1201) the constant a redefines the cosmo logical constant, and we assume that a = 0. The 
second parameter, fl, can be associated with the initial value of the Hubble function H(t ) = H : 
if we put Kf(H) from (121 into (125]), we obtain at t = t , that p=-3(l+aH 2 )+^. Then 
transforms into the bi-quadratic equation 



H — H \ Hn — 



kW \ kW 



3a 



3am 



a(t) 
a(t ) 



-37 







(29) 



which gives the following solutions 



nW s 
3aHl / 



± 



\4 



HI 



kW s 
3aHl / 



+ 



kW 
~3a~ 



a(t) 
a(t ) 



(30) 



Now we have to specify the signs of the parameters (3 and a. We assume here that a is positive and 



(3 is negative, thus, the function K,f(H) has one local maximum at H { 



(max) 



and two minimums 



6 



at Hi, 



± 



#-. We assume also, that the initial energy density W is restricted by inequality 



! (min) --y 2a 

kWo < 3cri/Q. Then the sign minus in front of the square root corresponds to imaginary solution for 
H(t). At a — > oo the second (real) solution yields the asymptotic value of the Hubble function 



Hi 



Hi 



kW 



3aH^ 



2 > 



Hl>\nl 



(31) 



When H = (^^ tt ) 4 , this parameter vanishes, i.e., H^ = 0. 

Now the parameter of acceleration can be presented in the form 



(4-3 7 )tf 2 (t) + (3 7 -2)tf ( 



2 

oc 



2[2H\t)-Hl 



(32) 



Clearly, the quantity —q(t) is non- negative for the traditional matter with 1 < 7 < ~, i.e., we deal 
with the accelerated expansion. When H^ 7^ 0, the acceleration parameter tends to one (— q(oo) = 1) 
in the asymptotic limit t — > 00, and the scale factor is of the de Sitter form a(t) oc explH^t]. When 
Hqq = 0, the solution for the Hubble function is 



H = Hn 



a(t) 
a(t ) 



(33) 



the scale factor a(t) takes the form 



a(t) = a(to) 



1 + - 7 ^o(t 



_4_ 

3i 



(34) 



and the acceleration parameter becomes constant 



_g(t)=l_- 7 , 

being, of course, non-negative for the traditional matter with 1 < 7 < ~. 
The effective (inertial) energy density u(t) (IT7|) takes now the form 



(35) 



u{t) = 3H 2 (t) - kW 



a(t) 
a(t ) 



-37 



(36) 



and tends asymptotically to the non- negative quantity uj (00) 
7r(t) (|T8l) satisfies the relation 



SH^. The effective (inertial) pressure 



ir{t)+u{t) 



-H\t)[l + q(t)]f\H) 



(37) 



In the asymptotic regime — q — > 1, thus, 7r(oo) = — cj(oo) = — SH^ is negative. In other words, the 
dynamic self-interaction influences the gravitating system analogously to the effective (depending on 
time) lambda-term, and this explicit analogy can explain the fact of accelerated expansion of the 
Universe. When the Universe reaches the de Sitter stage, this analogy becomes complete, since the 
constant value SH^ can be considered as effective (asymptotic) cosmological constant Aoo. 
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4 Conclusions 



1. We established a new model, which accounts for dynamic (inertial) self-interaction of gravitating 
systems. The extension of the theory of gravity is based on the introduction of some unknown 
function F(\l/ 2 ) depending on square of covariant derivative = ViUk of the velocity four-vector 
U h of the system as a whole. This unknown function should be determined using astrophysical 
and cosmological observations, as well as gravitational experiments. Since F(^/ 2 ) is the function of 
velocity derivatives, the Newtonian limit of this modified theory of gravity is not violated, i.e., all 
the supplementary terms in the Lagrangian decomposition start with (v/c) n , n > 1 and vanish at 

C — > CO. 

2. In this Letter we discuss only one example of exact solution to the extended master equations 
attributed to the isotropic cosmological model with vanishing cosmological constant. In the case, 
when the function F{^ 2 ) is quadratic in its argument, we have shown explicitly that the accelerated 
expansion of the Universe can be driven by matter with positive pressure (e.g., dust, ultrarelativistic 
fluid), when the dynamic (inertial) self-interaction of the gravitating system is taken into account. 

3. The dynamic self-interaction of the gravitating systems results in the appearance of additional 
effective (inertial) contributions to the total energy density u and pressure n (see (JT7I) . (|T8|) . (136]) and 
( |37|) ). These contributions look like the consequences of an effective lambda-term, which generally 
depends on time, but becomes constant on the de Sitter stage of the Universe evolution. In the 
asymptotic regime u becomes positive and ir negative, so that their sum vanishes. Thus, we can 
consider the accelerated expansion of the Universe to be driven by dynamic (inertial) self-interaction, 
which is characterized by negative pressure induced by back-reaction of the gravity field on the 
irregular motion of the system whole. 

4. In the context of isotropic cosmology only one type of the irregularity of motion, namely, expansion 
(VkU k 7^ 0), contributes into the dynamic self- interaction. For the static spherically symmetric 
gravitating systems the non-vanishing acceleration four-vector DUk should be taken into account. 
For the systems with pp-wave symmetry the shear tensor <Jik plays the main role in such dynamic 
self-interaction. As for rotating systems, such as spiral galaxies and pulsars, the contribution of the 
rotation tensor io^ seems to be the main one. 

5. The theory of dynamic self- interaction of an irregularly moving gravitating system should be 
experimentally tested. We could suggest at least three ways, how to do it. First, we intend to 
construct a realistic model of the Universe accelerated expansion by fitting the function —q[f(H(t))] 
(see (|27p ). based on the cosmological data; the reconstruction of the function f(H) will give us the 
basic function F(^ 2 ). Second, we hope to analyse post-Newtonian effects in this theory and to 
consider the corresponding experimental tests in the Solar system. Third, one can consider the flat 
rotation velocity curves of the spiral galaxies in the context of theory of dynamic self-interaction. 
We hope to study the mentioned models in the nearest future. 
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